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Abstract. A system of effective Einstein equations for spatially averaged scalar 
variables of inhomogeneous cosmological models can be solved by providing a 'cosmic 
CNJ ■ equation of state'. Recent efforts to explain Dark Energy focus on 'backreaction effects' 

of inhomogeneities on the effective evolution of cosmological parameters in our Hubble 
volume, avoiding a cosmological constant in the equation of state. In this Letter 
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C^) . it is argued that, if kincmatical backreaction effects are indeed of the order of the 

averaged density (or larger as needed for an accelerating domain of the Universe) , then 
■ the state of our regional Hubble volume would have to be in the vicinity of a far— 

, from-cquilibrium state that balances kinematical backreaction and average density. 
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This property, if interpreted globally, is shared by a stationary cosmos with effective 
O -1 equation of state p off = — 1/3 g cS . It is concluded that a confirmed explanation of Dark 

| | Energy by kinematical backreaction may imply a paradigmatic change of cosmology. 
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PACS numbers: 04.20.-q, 04.20.-Cv, 04.40.-b, 95.30.-k, 98.80.-Es, 98.80.-Jk 

1. Effective Einstein equations and the cosmic equation of state 

To set notation and to provide the framework for our argument, we recall a set of 
effective Einstein equations [2j. The argument presented can be carried over to studies 
of inhomogeneous cosmologies covering the Early Universe and radiation-dominated 
epochs with the help of the more general effective equations developed in 

1.1. Averaged equations 

For the sake of transparency we restrict ourselves to the matter model irrotational dust. 
Adopting a foliation of spacetime into flow-orthogonal hypersurfaces (which is possible 
for irrotational dust) with the 3-metric in the line-element ds 2 = —dt 2 + gij dX % dX\ 
we define spatial averaging of a scalar field f on a domain T> with volume Vp by: 

/ Jd'X *(t,X') , (1) 
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with J := ^det(gij), where is the metric of the spatial hyper surf aces, and X 1 are 
coordinates that are constant along flow lines. Following [2 we define an effective 
scale factor through the volume of a simply-connected domain T> in a hypersurface, 
normalized by the volume of the initial domain T>\, a© := (Vp/Vp.)^ 3 . Employing the 
fact that, for a restmass preserving domain D, volume averaging of a scalar function \I/ 
does not commute with its time-evolution, (d t ^) v — d t (^) v = (0) v — {^0) v , 

we can derive an effective equation for the spatially averaged expansion (8) v = ¥2. = 
3— = : 3Hx> (with an effective Hubble functional Hx> (an overdot denotes partial time- 
derivative): setting ^ = 9, inserting Raychaudhuri's evolution equation for d t 9 into the 
commutation rule above, and using the effective scale-factor a© one obtains: 

^+4kG(q) v -A=Q v ; (g) v = -^. (2) 
ax> a-p 

The first integral of (J2J) is directly given by averaging the Hamiltonian constraint: 

where the total restmass Mp, the averaged density, the averaged spatial Ricci scalar 
(72.) v and the kinematical backreaction term Qx> are domain-dependent and, except the 
mass, time-dependent functions. The backreaction source term is given by 

Q V : = 2 (II) V - 2 - (1)1 = \{{0- (6) v f) v - 2 (or 2 ) v ; (4) 

here, I = Q\ and II — |[ (O^) 2 — 0^0^] denote the principal scalar invariants of the 
expansion tensor, defined as minus the extrinsic curvature tensor, —Kij : = Qij. In the 
second equality above it was split into kinematical invariants through the decomposition 
Qij = \g-ij6 + o~ij, with the rate of expansion 9 = Q\ and the rate of shear a 2 = \oija l K 
(Note that vorticity is absent in the present gauge; we adopt the summation convention.) 

The time-derivative of the averaged Hamiltonian constraint Q agrees with the 
averaged Raychaudhuri equation (J2J) by virtue of the following integrability condition: 

d t Q v + GHvQv + d t (K) v + 2H V (K) v = & d t ( Q v a% ) + 4 d t ( (K) v a% ) = .(5) 
1.2. The cosmic quartet 

We may further introduce dimensionless average characteristics as follows j2]: 



Um - ml ' K - ml ' ml ' — ml ' [b) 

where we have employed the effective Hubble-functional Hx> introduced above that 
reduces to Hubble's function in the homogeneous-isotropic case. With these definitions 
the averaged Hamiltonian constraint (JHJ) reads: 

fi£ + nj + n£ + = 1 , (7) 

providing a scale-dependent cosmic quartet relating all relevant "cosmological 
parameters" . For Qv = the above functionals reduce to the corresponding parameters 
of the standard homogeneous-isotropic models. 
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The effective cosmological parameters defined in (JBJ) can be considered to provide a 
fair representation of the values which also an observer would measure in a sufficiently 
shallow survey region T> (the light-cone effect is not taken into account). We may 
therefore discuss estimates of those parameters in comparison with observed values. 
Note, however, that the interpretation of observations is mostly done by employing 
a standard Friedmannian cosmology as a 'fitting model' and therefore, geometrical 
inhomogeneities (that are hidden in the definition of the spatial averages in the 
Riemannian volume element, cf. Eq. are ignored [H]. 

1.3. The cosmic equation of state 

The above equations can formally be recast into standard zero-curvature Friedmann 
equations with new effective sources ([3]: Corollary 2)\: 

<&-^) v ~Qv~(n )v ; P?, = ~& + i ^a ■ (8) 

3— = A - 47rG(eS + 3P?s) ; 3fl£ = A + SiGgf, ; gf s + 3H V (<& + pj) = . (9) 

Eqs. © correspond to the equations (J2J), © and ©, respectively. In these equations 
the kinematical backreaction term Qd itself obeys a stiff equation of state mimicking a 
dilatonic field in the fluid analogy (for further implications see [3]). 

Given an equation of state in the form = /3(f?^,a£>) that relates the effective 
sources (jHJ), the effective Friedmann equations (jHJ) can be solved (one of the equations 
(0) is then redundant). Therefore, any question posed that is related to the evolution 
of scalar characteristics of inhomogeneous universe models may be "reduced" to finding 
the cosmic state on a given spatial scale. Although formally similar to the situation in 
Friedmannian cosmology, here the equation of state is dynamical and depends on details 
of the evolution of inhomogeneities. In general it describes non-equilibrium states. 



2. Explaining Dark Energy through kinematical backreaction 

The 'coincidence' that a Dark Energy source (modeled in the simplest case by a 
cosmological constant) starts to dominate around the epoch when also structure enters 
the non-linear regime suggests that there could be a physical relation between the 
effect of structure on the average expansion (known as backreaction effect) and the 
Dark Energy gap found in the Friedmannian standard model, thus providing a natural 
solution to this coincidence problem (besides the fact that more "exotic" explanations 
are then not needed). 

The averaging problem in cosmology has a long history including calculations of 
the backreaction effect shortly after George Ellis jH] has pointed out its importance ( [TT] 

X Note that in this representation of the effective equations p e s just denotes a formal "pressure" : in 
the perfect fluid case with an inhomogeneous pressure function the foliation has to be differently chosen 
and there is a further averaged pressure gradient term [3]. 



A cosmic equation of state for the inhomogeneous Universe 



4 



and many works thereafter; references may be found in [7HT31III1IH1EI)- The new input 
into this discussion is the (not unsupported) claim by Edward Kolb et al. ^3] that 
the condition (fTU|) below could be satisfied within our regional Hubble volume, hence 
providing a smart explanation of the Dark Energy problem ^3] without employing a 
cosmo logical constant, quintessence or corrections to Einstein's laws of gravity. 

Setting A = in Eq. the condition for an accelerating patch T> of the Universe 
directly follows: 

Q v > 4ttG (q) v ; (q) v = —\ > . (10) 

(With regard to (HJ), in order for Q v to be positive, expansion fluctuations would 
have to dominate over shear fluctuations.) This regional condition is weaker than 
the requirement of global acceleration, since it accounts for the regional nature of our 
observations. There is, however, a large body of opponents including myself who do 
not think that the condition (jlOj) can be met within the standard picture of structure 
formation from CDM initial conditions. A number of caveats would have to be overcome 
related to explicit calculations of kinematical backreaction and observational constraints, 
which both will be discussed in more detail in a forthcoming paper ^j. However, those 
caveats largely depend on assumptions that would consider perturbation theory on 
a Friedmannian background, would extrapolate ([TUjl to the global scale according to 
the cosmological principle, and also would ignore the difficulties in relating the model 
parameters (JBJ) to observations [9 . It is important to keep this disclaimer in mind in 
what follows. Restricting our attention to the universe model, as we shall do, there is 
no such caveat in our line of arguments. 

Let us give an example that illustrates how strong the condition ((TQJ) appears, if 
we straightly compare the model parameters © with current observations. We first 
rewrite the condition (fTTH) in terms of the dimensionless characteristics ©§: 

n v 

-n? > x ■ (ii) 

We have to be aware that, if (jllj) holds on some large domain V, which we may take to 
be as large as our observable Universe |13j . then Hamilton's constraint in the form (J7|) 
also implies: 

nj + nS > i - 2n£ , (12) 

showing that, for a low density parameter, we would need a substantial amount of 
negative curvature (positive f^) in the inhomogeneous model (not in the 'fitting model') 
on the domain T>, if we put the cosmological constant equal to zero. (To reconcile a 
small curvature parameter in this condition with a non-vanishing cosmological constant 
would need an even larger value of than that suggested by the 'concordance model' 
of about 0.7.) The fact that a large value of kinematical backreaction goes along with a 
substantial amount of average Ricci curvature has also been stressed and discussed by 

§ Note that, for a positive Qt>, fig is, by definition, negative. 
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Syksy Rasiinen The condition (|12j) . taken at face value (there are arguments why 
this could be naive), would contradict the widely agreed expectation that the curvature 
should be very small (the 'concordance model' assumes an exactly zero scalar curvature). 

We, here, do not attempt to respectively verify or falsify the condition ((TOJ) , but 
instead follow a line of arguments that assumes ()10|) to hold. That is to say, i/kinematical 
backreaction can indeed account for the Dark Energy gap found in Friedmannian 
cosmology, then we have to closely examine implications of the above condition. 

Let us first remark that the physical contents of (fTUj) implies a strongly fluctuating 
cosmos, roughly speaking: fluctuations encoded in Qx> have to be of the same order as 
the average density. This can only happen, if there is a strong coupling of kinematical 
backreaction to the averaged scalar curvature, even on the Hubble scale, for if this 
coupling is absent, © admits the special solution (11) v oc a^ 2 and Qx> oc a^ 6 ([2] 
App.B), i.e. averaged scalar curvature behaves as in a constant-curvature Friedmannian 
model with effective scale factor a©, and fluctuations decay with the square of the 
(expanding) volume capturing what we may call cosmic variance. This particular 
solution mirrors what we would expect from standard cosmology. It implies that the 
averaged density AtcG (g) v oc a^ 3 would substantially dominate over Q v oc a^, 6 in 
a globally expanding universe model after some time, even if we would start with the 
condition (JTHJ). Therefore, a strong coupling of Qt> to (Tl) v that changes the dependence 
on the effective scale factor sufficiently, is key to the explanation of Dark Energy through 
kinematical backreaction. 

Notwithstanding, if we assume such a strong coupling exists in a realistic universe 
model, and if we suppose that our Hubble volume accelerates due to the fact that Qv 
dominates over 4tiG (g) v , then we are entitled to say that a typical Hubble volume 
would correspond to a non-perturbative state in the vicinity of Qx> ~ 47rG (g) v , i.e. it 
would not correspond to a perturbative state in the vicinity of a model with Qx> ~ 0, as 
expected in the standard picture of small perturbations of a Friedmannian background. 

We are now going to identify this state in the effective equations. For this end let 
us now extend the spatial domain T> to the whole Riemannian manifold E, which we 
assume to be compact. The cosmological principle would extrapolate the condition |T0|) 
to the global scale. However, since we are assuming a strongly fluctuating cosmos, it is 
more appropriate to allow for other Hubble volumes that are slightly decelarating with 
Qt> < 4:7iG (g) v , so that a typical Hubble volume would reside in a state close to the 
balance condition Qx> = 47rG (£?}-£,||- This balance condition furnishes an example for 
such a state. Extrapolating this condition to the global scale implies with (J2J) that the 
global effective acceleration vanishes, and together with (jSJ) we face the global stationarity 
conditions: 

r 

g s = AixG(g) J] ; (1Z) E = 127rG(0> s -6i?| ; H s = — , C = const. , (13) 

as 

with the global kinematical backreaction Qs, the globally averaged 3-Ricci curvature 

|| In [5] a conservative estimate, based on current observational parameters, shows that such a cosmos 
provides room for at least 50 Hubble volumes. 
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{TV) j,, and the total restmass of the compact universe model M s . The second condition 
above has been obtained by eliminating the backreaction term in Q using the first 
condition. Eliminating instead the density source we can determine the constant C by 
evaluating ® at the initial time: C 2 := — \ The cosmic equation of 

state for a stationary cosmos can be obtained by inserting (fT3|) into (jEJ). We find^f: 

Peff = -3 QeS ■ (14) 



Now, taking the time-derivative of the two conditions in ()T3|) and employing restmass 
conservation on £, d t + 3H% {g) E = 0, we obtain evolution equations for the global 
kinematical backreaction and the globally averaged Ricci curvature. These evolution 
equations are solved by: 

Q s = — — ; = 3 + _ 3 > ( 15 ) 

which indeed points to a strong coupling between kinematical backreaction and averaged 
scalar curvature: the rate of decay of Qs is in proportion to and can therefore be of 
the same order as 47rG (o) s today. The total kinematical backreaction QyYy. = 4:7rGM-£ 
is a conserved quantity. Actually, the solution (|T5|) provides the first example of an exact 
solution of the effective Einstein equations with a non-trivial coupling of averaged scalar 
curvature to kinematical backreaction; (fT3j) solves the integrability condition (JSJ). 

As in Friedmannian cosmology the fate of the effectively stationary inhomogeneous 
cosmos is determined by initial conditions, for the stationary effective scale factor 
fl£ = o>s + C(t — U) can be restricted to a 'Big Bang model' by setting as = 0, or 
it can emerge ^0] from the effectively static cosmos as = const., which is a subcase of 
the stationary one by setting H s = ; C = <^> (^)s = 3Qs(tj). The globally 
static inhomogeneous cosmos is characterized by the cosmic equation of state: 

(ft) E = 3Q S = const. =4- Pe E ff = f? e E ff = . (16) 



3. Concluding remarks 

We argued that an explanation of Dark Energy by kinematical backreaction effects and, 
hence, requiring the condition (fTTH) to hold in our Hubble volume, implies a cosmos that 
is dominated by strong expansion fluctuations. We may speak of a far- from- equilibrium 
cosmic state in contrast to a perturbed Friedmannian state. We further argued that a 
system featuring strong fluctuations would conceivably, with some probability, create 
regional Hubble volumes in the vicinity of the stationary state ([15)1 . We have noted 
that, even if the Universe started in the vicinity of such a state, then it is important 
for the survival of the condition (jlOj) that the evolution implies strong coupling between 
kinematical backreaction and averaged scalar curvature. We have shown that the 
globally stationary state indeed conserves strong fluctuations: the key-property is that 

% It is interesting to compare this condition with the investigation of backreaction in inhomogeneous 
cosmon fields by Christof Wetterich JH1> m particular with the 'cosmon equation of state'. 
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<2s oc a^ 3 oc in large contrast to the case of a perturbed Friedmannian cosmos that 
would likely evolve into small fluctuations on the Hubble scale. 

The picture that emerges entails a regionally fluctuating cosmos, where the region 
is as large as our Hubble volume. In fact this would imply a paradigmatic change of 
cosmology, since the properties seen on the Hubble scale are not extrapolated to the 
global scale as in Friedmannian cosmology. This remark is more evident, if a globally 
static model is envisaged: a globally static, but inhomogeneous cosmos is conceivable 
without employing a compensating cosmological constant in contrast to the classical 
Einstein cosmos. This mathematical possibility (put into perspective in jl]) would attain 
the status of a viable physical model, if the "classical" explanation of the Dark Energy 
problem in terms of kinematical backreaction effects were true. In such a cosmos the 
averaged scalar curvature is, for a non-empty Universe, positive; the global kinematical 
backreaction term takes the role of a positive cosmological constant. Regionally, such 
a cosmos features exponential gravitational instabilities, i.e. it has strong matter and 
curvature fluctuations. 

In a forthcoming work [S] further implications are investigated. In particular, 
it is argued that, e.g. a globally static dust cosmos, provided it originates in this 
state, could be stabilized by backreaction and averaged curvature in contrast to the 
global stability properties of the classical Einstein cosmos £Q. It is further clarified 
that a fluctuation-dominated cosmos, expressed in a thermodynamic language, is in 
a far- from- equilibrium state compared with the Friedmannian "equilibrium state" by 
employing an entropy measure proposed in [12] . which vanishes for Friedmannian 
cosmologies ("zero structure") and is positive and time-dependent for the globally 
stationary, but inhomogeneous cosmos. 

This work was supported by the Sonderforschungsbereich SFB 375 'Astroparticle physics ' by the 
German science foundation DFG. Special thanks go to Stephane Colombi and Sabino Matarrese for 
stimulating discussions. 
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